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This  report  describes  the  achievements  of  the  research  program 
sponsored  by  AFOSR  under  Grant  AFOSR-80 -0029  during  the  period 
1  October  1979  to  1  October  1980.  The  object  of  the  research  program 
is  to  investigate  laser  interaction  with  dense  metallic  plasmas  and  its 
optical  as  well  as  X-ray  emissions.  We  study  the  dynamics  of  a  high 
density,  relatively  low  temperature  (20  eV)  plasma  when  illuminated 
by  a  high  intensity  CO^  laser  beam.  The  investigation  is  carried  out 
both  theoretically  and  by  means  of  an  advanced  numerical  simulation 
code.  The  following  sections  give  a  description  of  the  various  efforts 
and  accomplishment. 


A.  THERMAL  SELF -FOCUSING  OF  A  HIGH  INTENSITY  LASER  BEAM 
IN  A  HIGH  DENSITY,  LOW  TEMPERATURE  PLASMA 

Our  analysis  begins  by  constructing  the  governing  equation  for 

the  propagation  of  a  pulse  of  electromagnetic  wave  through  a  plasma. 

We  start  with  the  Maxwell  equation  in  a  general  dielectric  medium 

(scalar  dielectric  constant  e) 


(A  - 1) 


cVe  -  4^  +  Sr  * 

31  at2 


where  E  is  the  electric  field  in  the  radial  direction  and  j  is  the  corres¬ 
ponding  current  density.  By  means  of  the  transformation 
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we  remove  the  last  space  and  time  dependence  and  obtain  the  following 
equation  for  the  wave  envelope  4* : 


(A-2)  2iff  +  2i  =  ^  +  +f*- 

Equation  (A -2)  has  been  written  in  non-dimensional  coordinates  (kx,  eQi|h) 
the  nonlinear  term  f  is  given  by 


(A-2a) 
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n_(k  I  .  t) 

(«  -  — 5 - ) 
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Equation  (A-2)  can  be  put  into  standard  form  by  means  of  a  transformation 
to  wave  frame  coordinates  (?,  t),  where 


§  =  z  -  t. 

W e  then  obtain 

(A-3)  2i  vj;  +  ftp  . 

In  Eq.  (A-3)  the  nonlinear  term  f  is  both  space  and  time  dependent. 

We  apply  the  "moment"  method  to  the  analysis  of  Eq.  (A-3).  Its 
success  depends,  in  part,  upon  identifying  at  least  two  invariants  asso¬ 
ciated  with  Eq.  (A-3).  The  first  constant  of  the  motion  is 

(A-4)  I  =  J dx  dy 

The  second  constant  of  the  motion  is 

t  2 

(A-5)  I2  =  /  dx  dy  {  |Vi+  |2  -  /  dt'f2-^-}  . 

Using  these  constants  of  the  motion  we  can  construct  from  Eq.  (A-3) 
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a  relation  governing  the  behavior  of  the  pulse  radius.  First  we  define 
the  radius  based  upon  the  second  moment  of  the  radial  distribution 


(A -6)  <  a2>  =  J dx  dy  r^  |4>  | 2  /  J  dx  dy  |i|j  j  2. 

The  governing  relation  is  then  found  to  be 


(A-7) 


2  2 
d  <  a  > 


dt 


-  j  dxdy  { 1+  I2  + f l4»  |2  -  j  dt  f  }. 


In  order  to  progress  we  need  to  assume  a  radial  shape  function  for  the 
pulse  intensity: 

2  2*  2«U)  2  2 

k  I  =  — 2 - exp  /a 

a  (|.t) 

where  s(£  )  is  a  parameter.  Then  we  can  write  Eq.  (A-7)  in  the  form 

(A-8)  (a_  }  =  2a  2  -  4  J  (t,  a)  -8  J  dt'  J(t',a')— , 

at  u  a 

where  J  is  a  convolution  integral 

* 

00  2  2  2 
J(t,  a)  =  J  dr  f(r,  t)(l  -  “ )r  exp(-  r  /a  ). 

0  a 

Since  f  >  0  at  r  =  0,  and  f  —  0  as  r  —  oo,  we  can  show  that  J  is  always 
positive.  Thus,  as  the  pulse  begins  to  penetrate  the  plasma,  it  interacts 
with  it  and  creates  a  density  depression  (f  >0).  At  early  times  the  pulse 
radius  is  determined  by  the  relative  magnitude  of  the  two  terms 

2a^  -  4J(t,  a).  If  the  first  term  dominates,  the  pulse  diffracts 

2  2  2  2  2  2 
(d  (a  )/dt  >  0),  otherwise  it  will  self-focus  (d  (a  )/dt  <  0). 
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The  general  behavior  as  indicated  by  Eq.  (A -8)  conforms  to  our 
expectations.  We  expect  the  front  edge  of  the  pulse  to  be  diffracted 
while  the  later  parts  (experiencing  the  nonlinearity  created  by  the 
front  part)  undergo  more  focusing. 

To  proceed  further  we  need  to  specify  the  nature  of  the  interaction 
which  causes  the  nonlinearity.  The  response  of  the  plasma  density  to 
the  electromagnetic  field  needs  to  be  analyzed. 

The  laser  field  can  influence  the  plasma  target  in  several  ways: 
through  the  ponderomotive  force  effect  and  through  heating  the  plasma. 

In  our  low  temperature  high  density  plasma  the  collision  frequency  is 
sufficiently  large  as  to  make  plasma  heating  the  dominant  mechanism. 

At  the  same  time,  the  collision  frequency  is  sufficiently  small  as  to 
make  electron  thermal  conduction  (rather  than  electron-ion  collisional 
energy  transfer)  the  dominant  energy  dissipation  mechanism.  Thus 
the  temperature  equilibrates  at  a  much  faster  time  scale  than  the 
hydrodynamic  flow  time.  Hence  the  plasma  interacts  with  the  plasma 
mainly  through  the  temperature  equilibration  process  described  by  the 
energy  balance  equation: 

2 

e  n  v  _ 

(A-9)  V±  •(KeVi  Te)  =  eE  •  j_  =  - $—*1.  |  . 

m  Co 

e 

We  linearize  this  equation  and  combine  it  with  the  continuity  and  momentum 
equation  for  a  one -fluid  plasma.  The  electron  density  is  then  found  to 


satisfy  the  equation: 
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(A-10) 
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We  approximate  the  time  dependence  in  a  reasonable  manner  such  that 
this  enables  us  to  consider  the  radial  variation  separately.  Such  an 
approach  yields 


(A -11) 


2„  2  "l 

Cs  Vi  n 
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k  wp  2v  .  t  -(t-t1  )/t 

=  — - j  dt,|+|2e  R 
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where  t  is  a  characteristic  relaxation  time  of  the  flow,  e.  g. , 

XV 

r_  ~  a  /c  .  The  convolution  integral  J  in  Eq.  (A-8)  is  first  written  in 
R  —  o/  s 

the  form  (integrate  by  parts  twice): 


2  2 

up  a  i  oo  2  ? 

j  _  — 2 — —  J  dr(V  n)r  exp(-r  /a  ). 

4e  u  o 
o 


Using  (A -II)  we  now  find 

4  4  •  „  ,  2  2  ,tIv 

a  up  v  x  t  y  a  s(£  ) 

j  -  - 2__£ -  /  dt'  — - - - — -  exp{  -  (t-t’)/r  }  . 

32t re  u  KeT  o  {a  (4  ,  t)  +a  (|',t')} 
o  eo 

Here  £,  '  =  z  -  t'  ,  and  the  integral  with  respect  to  the  time  must  be  per¬ 
formed  at  constant  z  (rather  than  at  constant  |).  The  final  form  of  the 
self -focusing  equation  then  becomes: 
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d2(a2) 

dt2 


=  2a 


-  a 
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s<£ 1 ) 

a2+a>  2 


exp  {(t-t1  )/tr}  dt' 


o 


..  s(l"l . 

(a1  2+a"  ) 


exp  {(f  -t")/rR} 


where 

4  2  2  4 

£2  =  cop  v  .4*  a  /8tt  e  w  KeT 

o  ei  o  o  o  eo 

and  £  denotes  integration  with  respect  to  time  with  z  kept  constant. 

This  last  equation  will  be  solved  by  numerical  techniques.  It 
will  then  yield  the  variation  of  the  radius  of  the  laser  pulse  as  it 
penetrates  through  the  plasma  target.  Thus,  knowing  the  properties 
of  the  target  plasma  we  can  predict  the  behavior  of  the  laser  pulse  as 
it  penetrates  through  the  plasma. 


B.  DYNAMICS  OF  THE  STRONGLY  TURBULENT  PLASMA  GENERATED 
BY  LASER-PLASMA  INTERACTIONS 

A  high  intensity  laser  beam  disturbs  a  plasma  in  such  a  way  that 
instabilities  appear.  These  unstable  waves  grow  to  form  very  intense, 
spatially  localized  electric  fields.  Consequently,  the  plasma  is  ren¬ 
dered  strongly  turbulent.  Moreover,  such  localized  nonlinear  structures 
can  interact  with  themselves  as  well  as  with  individual  electrons  and 
ions,  thus  leading  to  a  significant  modification  of  the  zero  order  dynamic 
as  well  as  thermodynamic  properties  of  the  plasma  itself  (e.  g. ,  density. 


temperature,  flow  velocity,  etc.) 


A  model  has  been  developed  to  study  the  large  amplitude  kinetic 
Alfven  wave  in  more  than  one  dimension.  The  model  incorporated  the 
effects  of  finite  electron  inertia  and  finite  ion  Larmor  radius.  No  am¬ 
plitude  expansion  is  used  and  the  results  can  be  extended  to  arbitrary 
amplitudes  consistent  with  the  assumption  of  low  |3  in  the  plasma. 

Solutions  to  the  model  equations  have  been  analyzed  in  a  cylindri¬ 
cal  geometry.  These  solutions  take  the  form  of  a  waveguide  elongated 
in  the  direction  of  the  ambient  magnetic  field.  The  azimuthal  magnetic 
field  is  described  by  a  radially  varying  envelope  containing  sinusoidal 
oscillations  along  the  ambient  field  line.  The  solution  is  localized  trans¬ 
verse  to  the  background  magnetic  field.  The  width  of  the  radial  profile 

is  given  as  a  function  of  the  parameter  K  (=  V  /V  ).  As  K  in- 

z  A  phcl  S  0  z 

creases,  the  width  of  the  envelope  decreases  and  the  solution  amplitude 
increases  correspondingly.  In  the  opposite  limit,  as  decreases, 
towards  unity,  the  envelope  spreads  and  the  soliton  amplitude  decreases 
until  equals  unity,  when  the  nonlinearity  goes  to  zero.  At  this  point 
the  solution  spreads  toward  infinity  and  the  localization  effect  vanishes. 
We  call  this  the  plasma  wave  limit.  Therefore,  the  parameter 
characterizes  the  nonlinearity.  As  this  nonlinearity  increases,  wave 
steepening  is  observed.  However  the  width  of  the  envelope  reaches 
an  asymptotic  limit  of  approximately  three  sonic  ion  Larmor  radius. 

This  limit  suggests  that  there  is  a  stabilization  effect  against  collapse. 
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A  stability  analysis  is  performed  for  the  two-dimensional  cylindrical 
solitary  wave  solution  found  above.  First,  perturbations  preserving 
the  symmetry  of  the  solution  are  considered.  A  method  following  the 
work  of  Lyapunov  is  developed.  The  advantage  of  this  method  is  that 
it  can  be  applied  to  problems  where  an  explicit  analytical  solution  for 
the  "equilibrium"  state  cannot  be  found.  A  criterion  for  the  stability 
is  obtained  and,  in  pj  inciple,  can  be  computed  to  arbitrary  accuracy. 

This  criterion  is  then  applied  to  our  waveguide  solution  and  it  is  shown 
that  it  is  stable  against  symmetry  preserving  perturbations. 

This  part  of  the  research  constitutes  the  Ph.  D.  dissertation  of 


Mr.  J.  Sheerin. 
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C.  MODELS  OF  RADIATION  TRANSPORT  IN  DENSE  PLASMAS 

We  have  developed  a  radiation  transport  model  for  line  radiation 
whichfits  into  a  one -dimensional,  cylindrically  symmetric  Lagrangian, 

MHD  computer  code  with  nonequilibrium  ionization  dynamics.  A  laser 
heated  aluminum  exploding  wire  plasma  was  considered. 

The  model  combines  escape  and  transmission  probability  theory 
with  Eddington  theory  nd  is  henceforth  referred  to  as  the  EPDM  model. 
Empirical  functions  were  developed  to  relate  cylindrical  shell  proba¬ 
bilities  to  slab  probabilities  by  use  of  two  sets  of  model  problems: 
emission  from  an  infinite  homogeneous  purely  absorbing  cylinder  due 
to  an  axial  line  source  and  a  uniform  volumetric  source.  The  EPDM 
model  was  compared  with  the  variable  Eddington  theory  on  a  set  of 
problems  typical  of  the  ones  observed  in  the  plasma  simulation.  The 
EPDM  model  showed  similar  trends  for  all  sets  of  problems.  In  the 
optically  thin  or  purely  absorbing  regimes  the  EPDM  model  predicted 
higher  absorption  and  less  emission  than  the  variable  Eddington  calcu¬ 
lation.  However,  as  the  optically  thick,  scattering  dominated  regime 
was  approached,  the  EPDM  model  showed  lower  absorption  and  higher 
emission  than  the  Eddington  method.  This  was  to  be  expected  since 
the  approximations  made  in  the  EPDM  formulation  invalidate  it  in  this 
regime. 

The  EPDM  model  was  then  incorporated  in  the  magneto -hydrodynamic 
code.  The  results  of  a  standard  run  were  compared  to  those  of  a 
transparent  plasma.  The  most  dramatic  effect  was  that  the  addition  of 


10 


self -absorption  decreases  the  radiation  emission  from  the  plasma 
(see  Figure  for  a  typical  case).  The  emission  decreased  anywhere 
from  a  small  fraction  for  an  optically  thin  plasma  to  as  much  as  two 
to  five  times  for  the  optically  thick  plasma.  There  was  also  some 
redistribution  of  radiant  energy  due  to  the  fact  that  some  of  the  ex¬ 
cited  levels  could  decay  by  several  paths.  There  were  also  some 
minor  effects  on  plasma  dynamics.  The  electron  temperatures,  plasma 
radius,  and  the  effective  charge  were  greater  with  self -absorption 
while  the  ion  temperatures  and  the  densities  were  less  than  their 
respective  values  for  the  transparent  case. 

The  EPDM  model  is  now  being  used  within  the  MHD  code  to 
analyze  the  results  of  X-ray  measurements  of  the  emission  from 
laser-irradiated  exploding  aluminum  wire  plasmas  being  made  in 
our  laser -plasma  interaction  laboratory. 

D.  KINETIC  MODELS  OF  ULTRA  HIGH  DENSITY  PLASMAS 

In  the  theory  of  neutral  systems  such  as  gaseous  liquids,  the 
hard  sphere  fluid  model  has  served  as  a  calculational  basis  for  re¬ 
search  into  the  nature  of  real  fluids.  The  classical  one  component 
plasma  (OCP)  is  the  simplest  mathematical  model  of  a  system  incor¬ 
porating  Coulomb  interactions  and  serves  the  same  purpose  for  real 
strongly  coupled  plasmas  that  the  hard  sphere  model  has  served  for 
real  fluids.  The  model  consists  of  N  point  nuclei  of  charge  Ze  and 
mass  m  embedded  in  a  charge -neutralizing  fluid  (NZ  points  of  charge 
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Figure: 


Total  Emission  vs.  Laser  Energy  Deposition 
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-e  and  mass  m  smeared  into  a  continuous  bath), 
e 

To  further  specify  the  OCP  and  relate  the  model  to  real  physical 
plasmas,  we  need  to  define  a  few  parameters.  For  an  ion  fluid  of 
density  n  =  N/V  (V  is  the  plasma  volume),  an  appropriate  scale  of 
length,  the  ion  sphere  radius  (or  mean  ion  separation  distance),  is 
given  by 


,  4tt  -1/3 

a  -  <Tnl  • 

Because  of  the  1/r  nature  of  the  Coulomb  potential,  the  equilibrium 
properties  of  the  OCP  with  temperature  T  (=T^=Te)  which  deviate  from 
ideal  gas  behavior  are  governed  by 


r  =  (Ze)2  p/a 

where  P  =  (k_T)  * ;  is  Boltzmann' s  constant. 

D  O 

r  is  a  measure  of  the  amount  of  correlational  (potential)  energy 
in  the  plasma  relative  to  the  amount  of  kinetic  energy  present.  With 
this  relation  we  can  see  that  the  system  will  be  strongly  coupled  if  T 
is  of  order  one  or  greater. 

The  one  component  plasma  (OCP)  is  a  benchmark  problem  in  the 
study  of  systems  interacting  via  long  range  potentials.  The  model  has 
been  the  subject  of  much  investigation,  especially  in  the  last  decade, 
not  only  in  regard  to  the  understanding  of  abstract  systems,  but  also 
with  respect  to  a  more  accurate  description  of  real  physical  objects. 

In  particular  white  dwarf  stars,  neutron  stars  and  heavy  planet  cores 
are  objects  in  which  matter  is  in  a  state  similar  to  that  typified  by 


the  OCP  model.  The  plasma  produced  in  inertial  confinement  fusion 
experiments  can  be  dense  and  "cool"  enough  to  be  considered  strongly 
coupled. 

The  OCP  model  should  be  useful  in  interpreting  diagnostic  data 
from  such  experiments. 

One  of  the  primary  tools  in  plasma  diagnostics  is  the  scattering 
of  laser  radiation  from  the  plasma.  The  direct  interpretation  of  the 
scattered  light  data  is  accomplished  via  the  Fourier  transformed 
density -density  time  correlation  function,  which  in  this  case  describes 
the  differential  cross-section  of  electromagnetic  wave  scattering  from 
charged  particles. 

There  are  various  methods  available  for  calculating  this  dynamic 
structure  factor.  However,  few  of  these  methods  are  suited  for  strongly 
coupled  plasmas.  The  generalization  of  the  equation  of  hydrodynamics 
to  include  Coulomb  interactions  in  the  system  assumes  that  the  non¬ 
hydrodynamic  plasma  modes  are  collisionally  dominated,  an  assump¬ 
tion  that  ignores  the  independence  of  plasma  oscillations  from  conserved 
variables.  Hence  a  microscopic  kinetic  description  of  the  full  phase 
space  density  is  required  to  develop  a  first  principle  theory  applicable 
to  a  plasma  in  a  strongly  coupled  state. 

The  well  known  perturbation  expansions  and  low  density  kinetic 
approaches  to  obtaining  a  kinetic  equation  for  the  phase  space  density 
are  inadequate  for  strongly  coupled  systems.  A  somewhat  better  suited 
approach  involves  transforming  the  Liouville  equation  into  a  form  of  a 
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Figure 


Comparison  of  the  TCP  model  with 
physical  interest. 
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generalized  Langevin  equation  that  provides  a  point  of  departure  for  fur¬ 
ther  approximation.  The  result  is  a  kinetic  equation  involving  a  mem¬ 
ory  function  or  kernel  containing  the  static  correlation  information 
as  well  as  the  difficult  collisional  dynamics. 

The  separation  of  the  static  and  collisional  parts  of  the  memory 
function  and  the  ability  to  choose  a  collision  term  that  dynamically 
evolves  from  the  exact  static  correlations  statics  allow  one  to  utilize 
the  recently  available  "exact"  Monte  Carlo  (MC)  equilibrium  data  direc¬ 
tly  in  the  solution.  In  the  cases  where  the  computer  data  was  compared 
with  a  Debye -Hvlckel  equilibrium  correlation,  exact  in  the  long  wave¬ 
length  and  small  coupling  regimes,  the  consequent  dynamic  structure 
factor  was  found  to  be  in  closer  agreement  with  "exact"  particle  dynamics 
simulation  results  when  the  MC  correlations  were  used.  In  view  of  the 
short  wavelength  difference  between  the  Debye-Htickel  and  MC  tabula¬ 
tions,  this  conclusion  should  be  expected.  But  the  point  is  that  since 
equilibrium  data  is  fairly  easy  to  generate  (as  opposed  to  wavelength- 
frequency -coupling  specific  molecular  dynamics  simulations),  a  kinetic 
representation  should  make  full  use  of  such  information. 

A  variety  of  collision  models  have  been  developed  within  the  ki¬ 
netic  equation  framework  and  applied  to  the  analysis  of  the  OCP.  All 
are  capable  of  using  the  computer -evaluated  static  correlations,  and 
all  allow  the  kinetic  equation  to  be  solved  exactly.  In  the  course  of 


these  solutions  we  found  that  the  mean  field  or  Vlasov  approximation, 
which  completely  ignores  the  dynamics  of  the  memory  function,  is  not 
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adequate  for  long  wavelengths  or  large  coupling  (r).  The  Fokker  -Planck - 
like  Lenard-  Balescu  model  allows  the  inclusion  of  dynamics  in  the  sim¬ 
ple  form  of  a  collision  time.  As  long  as  small  wave  numbers  are  avoided 
and  the  collision  frequency  is  derived  from  available  molecular  dynamics 
data,  the  spectrum  of  S(k,  oj)  is  improved  dramatically.  The  Fokker- 
Planck  approximation,  however,  does  not  satisfy  the  collisional  invari¬ 
ant  of  momentum,  hence  the  anomalous  results  at  long  wavelength. 

The  exact  calculation  of  the  collision  term  to  second  order  in  fre¬ 
quency  modified  by  the  inclusion  of  two  unknown  relaxation  times  —  a 
model  that  does  satisfy  conservation  of  momentum  —  is  relatively  easy 
to  accomplish.  This  modified  Duderstar’t-Akcasu-Linnebur  (DAL)  model 
makes  use  of  more  known  computer  information  (via  the  "exact"  shear 
viscosity)  and  gives  results  that  are  in  very  good  agreement  with  simu¬ 
lation  calculations  over  a  wide  range  of  wave  number  and  coupling* 

The  Jhon -Forster  (JF)  model  is  an  extension  of  the  DAL  model 
to  include  the  energy  collision  invariant.  Due  to  the  fact  that  the  OCP 
is  dominated  by  the  plasma  oscillations  and  that  charge  and  mass  fluc¬ 
tuations  in  the  OCP  are  equivalent,  the  energy  modes  do  not  couple 
efficiently  to  the  mass  density  modes.  The  consequence  is  that  the  JF 
model  provides  only  slight  improvement  over  the  DAL  model  except  at 
very  large  coupling  and  very  small  wave  number.  The  amount  of  cal- 
culational  effort  expended  in  procuring  this  slight  advantage  along  with 
the  need  to  obtain  yet  another  transport  coefficient  and  still  unavailable 
static  correlations  indicate  that  the  JF  solution  is  more  extensively 
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approximated  and  not  as  efficient  as  the  DAL  calculation.  The  results 
can  even  be  worse  for  shorter  wavelengths. 

The  dynamical  collision  models  require  known  presumed  exact 
transport  coefficients  as  input  to  insure  correct  behavior  or  to  cover 
the  low  frequency,  low  wave  number  limit  of  an  interpolative  model. 
Thus,  no  information  on  these  coefficients  can  be  gained  by  utilizing 
these  terms  in  a  kinetic  approach  toward  a  calculation  of  transport 
properties.  In  addition,  for  the  strongly  coupled  OCP  we  found  that 
the  long  wavelength,  long  time  values  of  S(k,  w)  calculated  from  the 
four  collision  models  were  not  always  in  good  agreement  with  the 


temporal  simulation. 
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